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Dynamics of tunneling centers (TC) in metallic systems is 
studied, using the technique of bosonization. The interaction 
of the TC with the conduction electrons of the metal involves 
two processes, namely, the screening of the TC by electrons, 
and the so-called electron assisted tunneling. The presence of 
the latter process leads to a different form of the renormalized 
tunneling frequency of the TC, and the tunneling motion is 
damped with a temperature dependent relaxation rate. As 
the temperature is lowered, the relaxation rate per temper- 
ature shows a steep rise as opposed to that in the absence 
of electron assisted process. It is expected that this behav- 
ior should be observed at very low temperatures in a careful 
experiment. The present work thus tries to go beyond the ex- 
isting work on the dynamics of a two-level system in metals, 
by treating the electron assisted process. 

PACS indices : 66.30. Jt ; 72.10.Fk ; 05.30.-d 



I. INTRODUCTION 

The existence of double-well structures in metals has 
now been known for a long time. They arise in di- 
verse physical situations, the most well-known of which 
is metallic glasses. The existence of tunneling states in 
glasses was proposed long back [jy. It is believed that 
because the atoms in a glass are 'quenched' in a ran- 
dom configuration, there is a possibility of an atom, or 
group of atoms to see a complicated potential landscape 
in its neighbourhood, which may have at least two neigh- 
bouring local minima. In such a situation, the atom, 
or the group of atoms can spontaneously tunnel from 
one well to the other. In this, and similar situations, 
the entity is referred to as a tunneling center (TC). An- 
other well known, and well studied example is that of 
hydrogen trapped in Nb(OH) x . At very low temper- 
atures (~ 70 if), the hydrogen, which can normally dif- 
fuse quantum-mechanically, is localized in some trap-sites 
formed by oxygen atoms The geometry of the sys- 
tem is such that the hydrogen can stay in two close-by, 
energetically equivalent sites around the trapping atom. 
This results in the formation of a TC. Incoherent neutron 
scattering reveals the tunneling of the hydrogen, which 
is influenced by the conduction electrons of the metal. 
These experiments have been theoretically analyzed us- 
ing the model of a particle in a double-well, in contact 
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with free electrons [|| |H| . Lately, the conductance be- 
havior of narrow metal constrictions has been explained 
within the framework of scattering of conduction elec- 
trons from a TC ||. The behavior of certain heavily 
doped conducting polymers has also been attributed to 
the existence of TCs [Q. Universal conductance fluctu- 
ations have also been used to study the dynamics of a 
single two- level defect in a metal ||. In short, tunnel- 
ing centers in metals have proved to be of considerable 
experimental and theoretical importance. 

Theoretical scenario is the following. The TC can be 
classified to be either slow or fast, depending on its de- 
tailed behavior. When the TC is slow, the effect of the 
conduction electrons is to provide some kind of friction to 
the tunneling motion of the particle. What actually hap- 
pens is that the particle being charged, has a screening 
conduction electron cloud around it. As it tunnels from 
one site to the other, the electrons have to follow its mo- 
tion adiabatically. Now it turns out that the electrons 
cannot actually follow the motion of the particle instan- 
taneously - the adiabatic approximation breaks down. 
The tunneling of the atom is a local perturbation which 
can give rise to an infinite number of electron-hole pairs 
in the electron gas, with infinitisimal energy, and thus to 
the so-called infra-red divergence. The result is that the 
tunneling of the particle is governed by the relaxation of 
the conduction electron cloud around it. The tunneling 
is thus heavily suppressed and at zero temperature the 
particle tends to localize in a particular well. This small 
energy behavior of the electron gas can be described in 
terms of charge density excitations of the electron gas, 
which are approximately bosonic in nature. In addition, 
if the low-temperature dynamics of the particle is con- 
fined to the two-dimensional Hilbert space spanned by 
the two ground states of the isolated wells, one can map 
the problem onto the model of dissipative two-level sys- 
tem (TLS) studied extensively by Leggett et al ||. One 
has to choose what is called the Ohmic dissipation in 
order to correctly describe the behavior of electronic ex- 
citations. 

In a slow TC, the atom remains stationary during an 
electron scattering. On the other hand, if the TC is fast, 
the process where there is a scattering of an electron si- 
multaneously with the tunneling of the particle, has to 
be taken into account. This process was originally pro- 
posed by Kondo ]lO[ and further work along this direction 
has been done by Zawadowsky and collaborators [p~T|— p~5|| . 
Here the electrons assist the tunneling motion of the par- 
ticle, and hence the particle cannot be localized in one 
of the two wells, even at very low temperatures. The 
particle, tunneling from one well to the other, scatters 
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the electrons between different orbital angular momen- 
tum states. If we describe the TLS, formed by the parti- 
cle in a double-well, by a pseudo-spin 1/2, the situation 
is quite similar to the Kondo problem where the impu- 
rity spin-flip causes a spin-flip of the electrons. Here the 
role of the electron-spin in the Kondo problem is played 
by the different orbital angular momentum states of the 
electrons. In the present problem, the role of the electron 
spin itself is rather passive, and the two spin states act 
like two independent channels of electrons to which the 
TC couples. It has been shown earlier that this prob- 
lem can be mapped on to that of the two-channel Kondo 
model, originally proposed by Nozieres and Blandin [ |l6| , 
which has come under a lot of attention because of its 
non-Fermi-liquid low temperature behavior. 

In an actual TC in metals, both the processes should be 
present - infact, the amplitude of screening by conduction 
electrons is much larger than that for electron assisted 
tunneling (because the tunneling particle is heavy). The 
interplay between these two processes may yield inter- 
esting dynamics at low temperatures. The dynamics of 
a TLS, with electron screening, has been studied in de- 
tail, and the high as well as low temperature behavior is 
well understood. At very low temperatures, the TLS un- 
dergoes weakly-damped coherent oscillations, with a tun- 
nel frequency which is strongly renormalized by the cou- 
pling to the electrons. The damping coefficient increases 
roughly linearly with temperature, whereas the effective 
tunnel frequency follows a weak power-law. As the tem- 
perature is increased, the thermally excited electron-hole 
pairs destroy the coherence of the particle. As a result, 
tunneling gets incoherent, and resembles stochastic jump 
diffusion. The jump rate, from one site to the other, how- 
ever, has a very non-trivial temperature dependence. It 
follows a power-law, the jump rate decreasing with in- 
creasing T. 

The dynamics of a slow TLS, in all temperature 
regimes, has been studied within the so-called dilute 
bounce gas approximation (DBGA), to the underlying 
functional integral expression for the dynamics of the 
TLS §,11 . For weak coupling of the TLS to the elec- 
trons, some calculations have gone beyond the DBGA, 
to study the low-temperature dynamics HJ17J . These two 
approaches do not treat the electron assisted tunneling 
process. Very recently conformal field theory has been 
used to study the very low temperature scaling behavior 
of the two-channel Kondo model [Q. The results ob- 
tained are relevant for the anomalous resistivity of the 
metal arising out of the scattering of electrons off the 
TLS. This approach has been highly successful. We are, 
however, interested in the actual dynamics of a TLS and 
not on its effect on the metallic system. Such interest 
stems from the fact that while experimentally probing 
the dynamics of a TLS, one usually deals with things 
like the tunnel frequency and the relaxation rate of the 
tunneling atom. 

In this paper, we calculate the Laplace transform of 
the time correlation function of the TLS, which, in the 



context of neutron scattering from the tunneling parti- 
cle, is related to the incoherent structure factor. The 
correlation function will give the necessary information 
regarding the interplay between the tunneling dynam- 
ics and the relaxation behavior, so commonly seen in all 
quantum dissipative systems. 



II. THE MODEL 

In the following we assume that the TC consists of a 
single particle in an asymmetric double well potential. 
The temperature is low enough for the dynamics of the 
particle to be effectively described by the Hilbert space 
spanned by the two lowest states of the isolated wells. In 
this approximation, the TC becomes a TLS. Now, the 
interaction of the TLS with the conduction electrons can 
be described by the following Hamiltonian (see e.g., | ff2| ) 

H' = 7^eer 2 + iftA er x 
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In (|l|), TLS is describe by the pseudo-spin a. The two 
eigenstates of o~ z , | ± > correspond to the tunneling par- 
ticle being in the left and the right wells, respectively. 
The electrons are described by a spherical wave repre- 
sentation, due to the presence of a point scatterer in the 
Fermi sea. The first two terms in (Q) describe the bare 
TLS, e being the asymmetry between the two wells and 
Ao the bare tunneling frequency. The third term repre- 
sents the screening of the TC by the conduction electrons. 
The fourth term in ([[]) describes electron assisted tunnel- 
ing, Ai being the strength of the process. The last term 
is the Hamiltonian of the electrons in the partial wave 
representation. The indices k and k! represent the wave- 
vectors of the electrons; I is the index for the partial wave 
and in the present model runs over two values, + and — ; 
p represents the spin of the electrons and runs over two 
values ±1/2; e p (k) is the energy of the electrons with 
wave- vector k and spin p; p is the density of states at the 
Fermi level of the electron gas; K is a dimensionless cou- 
pling constant which describes the screening interaction 
of the TLS with the electrons, whose value is restricted 
to be < K < 0.5 Operator c pk+ (c pk+ ) creates 

(annihilates) an electron in a state with spin p, a wave- 
vector k, energy e p (k) and partial wave state +. The 
energy of the excitations is assumed to have an upper 
cutoff, equal to the Fermi energy of the electron gas, cf- 
The Hamiltonian in ([!]) can, infact, be thought of as a 
highly anisotropic two-channel Kondo Hamiltonian with 



2 



two external magnetic fields e and Ao, pointing along 
z- and x-axis respectively. It also differs crucially from 
the Kondo problem due the fact that the value of the 
so-called Jh, in the present case is very different. Due 
to this, the present model shows coherent oscillation of 
the pseudo-spin, which is absent in in the conventional 
model. 

We have shown earlier that the dynamics of the sin- 
gle channel Kondo problem can be efficiently studied by 
bosonizing the Hamiltonian ^0|. The success of this 
method in treating the effect of the low-energy excita- 
tions on the impurity spin motivates us to also apply it 
to the present problem. The two spin channels of the elec- 
trons being completely independent, the bosonization of 
one does not affect that of the other. Consequently, the 
bosonized form for ([!]) can be written down in a manner 
similar to that of the single channel Kondo Hamiltonian 
(for a detailed account of bosonization the reader is re- 
ferred to The resulting Hamiltonian then assumes 
the following form 

H' = -hea z + ^hA cr x 

+ lhy/2nKv F /L<r z J2 vW* / *''" ( b tk + b pk) 

k,p 

+ -HAi ^[tr+e - ** + T-e +ip ] + nuj kb f pk b pk (2) 

p kp 

where £ p = \ J2 P ( nv F / 'Luk) 1 / 2 exp(~uj k /2uj c )(b pk - b pk ). 
In (||), p is the channel index, cu k = vpk, vp being the 
Fermi velocity and tu c is a high frequency cutoff related 
to tp . L is the length of the normalization box and the 
limit L — > oo is taken such that (2n/L) J^k ~ ¥ I dk. The 
bose operator b L creates a particle in the state k of the 
p'th channel, which is an excitation of the electron gas. 
The charge density excitations drop out of the problem. 

If Ai is zero, the Hamiltonian given by (g), with- 
out Ao, can be diagonalized exactly. The strategy 
then, is to treat the terms proportional to A and Ai 
as a perturbation. To this end we perform a unitary 
transfor matio n UH'U^ 1 on the Hamiltonian where U = 
exp^ y/K/2 ^2 p £ p ] . The transformed Hamiltonian has 
the following form 

H = -%ea z + ^ hu k b pk b pk 

k,p 

+~hA Q [<j + e^^ + ^ + «7_e->/?K 1+ &>] + -TiAj 

x J2lv+^ iil+i2) ~ ip + <T_e-V?« 1+ ^)+^] (3) 
p 

One can see that for Ao, Ai = 0, the Hamiltonian in (||) 
is diagonal in the usual representation. In the following 
we denote the first and the second terms of (||) by Hg 
and Hb, respectively, and the other two terms by Hi. A 
suitable perturbation theory can now be done on H by 



treating Hj perturbatively, which amounts to treating K 
exactly and doing an expansion in Ao and Ai. 

The quantity of interest here is the Laplace trans- 
form of the time correlation function of the pseudo-spin, 
C(t) =< a z {0)a z {t) >, given by 

C(z) = / e~ zt < a z (Q)a z (t) > dt. (4) 

Here, the angular brackets denote canonical ensemble av- 
erage, and the Heisenberg time evolution of a z is dictated 
by H. In the following analysis, the Greek indices fi, v 
denote the impurity spin states, and the "states" of the 
Liouville operators are denoted by | v\l) etc. We shall 
now focus our attention on C(z) which can be put in the 
following form: 

C(z) = Y<v\p s \vXv\S z \v> 

I/, v' 

{y,v\ [U(z)] av | v',v') <v'\S z \v' > (5) 

In writing (^|) we have made use of the Liouville operator 
formalism to introduce a "bath-averaged time evolution 
operator" [U(t)] av , where U(t) = e lLt , L being the Li- 
ouvillian associated with H. In addition, we have fac- 
torized the canonical density matrix p, as p ~ Ps ' Pb, 
where ps is the density operator associated with H$ 
and ps is the density operator for a bath of noninter- 
acting bosons given by Hb- The bath average [...] au 
implies a multiplication by ps and a trace over the 
bath states. The bath-averaged time evolution opera- 
tor contains all the information regarding the influence 
of the electronic environment on the TLS. We calcu- 
late [£7(z)] au using the resolvent expansion formalism 
where the resolvent is treated perturbatively to yield a 
"self-energy" which is second order in Hi, resulting in 
[U{z)]av » [z - iL s + {Li{z - iL s - iL B )~ 1 L I } av }- 1 , 
where L§, Li and Lb are Liouville operators associated 
with Hs, Hi and Hb respectively gj. The calculation 
proceeds along the lines of Ref. || which treats the dy- 
namics of a TLS in contact with a bosonic bath. 

We first calculate the matrix for the self-energy [Li(z— 
iLs — iLs) 1 Li] av and it turns out that in the second or- 
der resolvent expansion, it can be represented in terms of 
certain correlation functions of a gas of free bosons. In- 
fact, one needs quantities like < e +a ^ 1 ■e~ b ^ 1 ^ >b where 
the subscript B indicates that the correlation function is 
calculated using only a free-boson Hamiltonian |3| . These 
correlation functions can be calculated using certain well 
known properties of a set of harmonic oscillators using 
the approximation T <C ?kJ c /fcs and t <C 1/oj c ||. The 
self-energy, being related to the Liouville operators, is a 
4x4 matrix within the space of the states of the pseudo- 
spin. But this matrix turns out to be block diagonal, and 
the 2x2 blocks can be handled with ease. Consequently, 
the block-diagonal matrix for the self-energy is combined 
with the diagonal matrix for (z — iLs) and inverted to 
yield the averaged time-evolution operator [U(z)] av . 
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III. RESULT 

Once the bath-averaged time-evolution operator is 
known, we are all set to calculate the Laplace trans- 
formed correlation function, the quantity of central im- 
portance to us. We plug in the form of [U (z)] av in (||) to 
yield : 

A ( s z + tanh (fie/3/2) + ie) - $_(z - ie)} 

[Z> z[z + $+(z + ie) + $+(z-ie)] ' 1 ' 

where $>±(z) — F(z) ± F'(z), and 

F( Z ) 1 a*(o/u, ^K m-mm + K+z/e) K 
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-z/Q) \ 


A Ai r(i + 25)r(l - <5 


+ z/Q) 



4tt 2 (z-Se)T(l + S + z/Q) \w, 
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(7) 



In (0), 6 = 2irk B T/h, r) = V2K - K, 6 = V2K - 2K 
and is the same as F(z) except that the exponential 
factors are replaced by their complex conjugates. Some 
comments at this stage. The expression for the Laplace 
transformed correlation function, given by (Q), describes 
the dynamics of a TC in a metal. The dynamics is gov- 
erned by the screening of the TC by the conduction elec- 
trons, as well as the electron assisted process. The first 
term represents the effect of screening of the TC by the 
conduction electrons, while the second term accounts for 
the electron assisted process. The third and the fourth 
terms in ^ are much smaller in magnitude compared to 
the first two, because of the presence of the extra factors 



IV. DISCUSSION AND CONCLUSION 

The second order treatment of Ao and Ai is equivalent 
to some kind of a DBGA if one were to employ the func- 
tional integral formalism, pioneered by Leggett et al 
Let us, for the time being, neglect electron assisted tun- 
neling, by putting Ai — 0, which reduces the expression 
in (pi) to the following 



C(z) 



where 



z + i tanh (^hefl) sin ttK[G(z + ie) — G(z — ie)} 



z[z + costtK{G(z + ie) + G(z - ie)}} 



(8) 



G{z) - -A (O/c c ) {z + eKni _ K + z/Q y (9) 



This is just the DBGA result for a dissipative asymmet- 
ric TLS, for ohmic dissipation, derived by Leggett et al 
B using the path-integral formalism, and by Dattagupta 
et al H using the resolvent expansion formalism. By in- 
cluding the Ai-dependent terms in the expression for the 
correlation function, we extend the DBGA result to the 
case where the electron assisted process plays an impor- 
tant role, in addition to the screening effect. 

In order to make explicit analysis simpler we look at 
the particular case of a symmetric double- well where e — 
0. The real part of C(z), which may describe the neutron 
scattering structure factor, is in general not a Lorentzian. 
In order to force a Lorentzian form, we further assume 
that K <S 1. With this simplification, the expression (||) 
can be approximated by 

C(Z) » =r- sr- , (10) 

2 + A 2 ) /(z + 70 ) + 2A 2 (z + 7l )' 

where A = A (9/w c ) K [r(l - 2K)Cos(irK)} 1 / 2 , Ai = 
{A 1 /Lu c )(io c /e)^r{2 V )Co S (nn)Ml - 2TJ)} 1 / 2 , 7o = K@ 
and 71 = 77O. In writing (HG) we have neglected the 
last two terms in (Q) (recall the approximation uj c 3> 
0). The real part of the Laplace transformed corre- 
lation function given by ( |To| ) describes two Lorentzian 
lines of width 7£ = 7o + 2A 2 7 i/(l + 2A 2 ), centered at 
w = ±[A 2 /(l + 2A 2 )- 7o + 2A 2 7l /(l + 2A 2 )] 2 . One can 
see that the tunnel frequency of the particle is renormal- 
ized due to the effect of metallic electrons. The renormal- 
ization is not only due to the overlap of electron states 
corresponding to the two positions of the particle (which 
leads to a reduction factor of the form {UbT /hio c ) K ) but 
also due to the effect of electron assisted process. This 
will show up in the effective tunnel splitting of a particle, 
at very low temperatures where it is in the coherent tun- 
neling regime. The width of the peaks, because of the Ai 
dependent term, also deviates from the usual "Korringa 
form" at very low temperatures. 

The DBGA is known to break down at very low tem- 
peratures for small frequencies. In order to circum- 
vent this difficulty we shall concentrate on the real part 
of C(z) near the resonances. Without any interac- 
tion with electrons the resonances for the TLS are at 
2 = ±iAo- So, we replace F(z) in the expression (0) 
for C(z), by F(iAo). Thus, the relaxation rate is given 
by 7 = Re[G(iAo)], and the tunnel-splitting is given by 
lo = —Im[C(iAa). Figures [i]and ^ show plots of the two 
against temperature. In the plots, everything is scaled 
with Ao, to yield dimensionless variables. One can see 
that for low enough temperature, the electron assisted 
process has a dramatic effect on the relaxation rate. At 
very low temperature, such a system is known to go into a 
Kondo-like correlated state. This may be the reason for 
this increased contribution to relaxation. This feature 
can possibly be experimentally observed. The tunnel- 
splitting, however, does not show a significant qualita- 
tive change in its temperature dependence. The plots 
indicate that at high temperatures the electron-assisted 



4 



process does not have a significant effect. But as the tem- 
perature is lowered, it is this process which dominates the 
dynamics of the TLS. 

In conclusion, we have studied the dynamics of a parti- 
cle in a double-well potential, in the presence of conduc- 
tion electrons where the effect of electron assisted tun- 
neling is taken into account in addition to conventional 
damping effects. In the appropriate limit the relevant 
expression reduces to the DBGA result for the dynamics 
of an Ohmic dissipative two-level system. In the coher- 
ent tunneling regime, the tunnel frequency is modified by 
the effect of electrons. The relaxation rate shows an un- 
expected rise as the temperature is lowered, as opposed 
to the case where the electron assisted process is absent. 
We believe that this feature should be observable in a 
careful experimental study. 
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